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SUMMARY 

Equations expressing the heat transfer and evaporation from wetted 
surfaces during ice prevention have been simplified and regrouped to 
permit solutions by simple graphical means. Working charts for quick 
and accurate anti-icing calculations are also included. 


INTRODUCTION 

Solution of the general problem of heat and mass transfer from a 
wetted surface in forced convection is quite involved and tedious. 

The calculations that often must be made point by point along a surface 
are complex even after the basic variables and ambient conditions have 
been determined. The method of solution customarily performed involves 
several trial-and-error calculations that are intermediate 'between the 
final answer and the basic factors that define a particular anti-icing 
situation. 

This investigation was conducted at the NACA Lewis laboratory to 
simplify the method of calculating heat and mass transfer after the basic 
factors (such as ambient temperature, rate of water interception, and 
relative velocity) are known. Calculations will be greatly simplified by 
a fundamental rearrangement of terms in the conventional equations, which 
will eliminate the trial-and-error calculations and permit a rapid 
graphical solution. The solution will be limited to the normal range of 
anti-icing conditions and to the case of surface temperatures above 
32° F with liquid-water interception by the surface. 


DEVELOPMENT OF GRAPHICAL SOLUTION 

Simplification of Conventional Equations 

The solutions most frequently desired in anti-icing calculations 
are for the heat-transfer rate associated with prescribed anti-icing 
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conditions and the corresponding evaporation rate. These terms consti- 
tute the two principal solutions of the anti-icing problem and are 
developed separately. 

Heat tr ans fer. - The heat transfer in forced convection from a sur- 
face subjected to water impingement and heated above 32° F has been rep- 
resented by the total heat transferred by convection and by evaporation 
of surface water, plus the sensible heat change of the impinging water. 
These heat-transfer processes are discussed in detail in references 1 to 
5 and are summarized in the following five conventional equations: 

q = h a (tg - t^) + hg^X - l) K (tg - t^) + My-Cy^tg - ty) (l) 

where 



e 0 “oT 0 

3.49 + 9.05 — +0.000683 — — 

\ P0 p 0 


' L i e z L o e o 
TqPo 


(The symbols used herein are given in the appendix. ) 

A solution of equation (l) requires the prior calculation of equa- 
tions (2) to (5) • Because of the dependence of e d on t d and of e^ 

on t-,, equations (2) and (5) are both trial-and-error solutions. In 
order to eliminate these laborious calculations, a single equation will 
be derived conta ining only a minimum number of basic variables. 

The n umb er of basic variables that define an anti-icing situation 
can be conveniently reduced to five: t s , tQ, Vq, Pq, and p-j ■ . In addi 

tion, three complex terms that cannot at the present time be -generally . 
reduced to basic variables must be retained in the final expression. 
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These terms hg, My., and K are discussed in a subsequent section of 

this report. In order to obtain a simplified relation using the pre- 
ceeding eight terms; the followin g assumptions are made: 

(1) In the air flow over a body, temperatures and pressures change 
adiabatically too quickly to permit any change in the state of entrained 
moisture; the local-stream vapor pressure (as shown in reference 2) is 
then given by 



( 6 ) 


(2) The local velocity and pressure may be related by the incom- 
pressible dry-air relation 

v i 2 - v o 2 - £- (pi - p 0 > 7 °- 7 (n 

Further simplification is obtained by expressing the following minor 
variables as constants near the middle of their expected range in anti- 
icing calculations: 

(1) Recovery factor R, assumed as 0-85, corresponding to the con- 

servative case of laminar flow 

(2) Latent heat of vaporization L, assumed as 1066 Btu per pound, 

corresponding to an average surface temperature of 50° F 

(3) Ambient absolute temperature (for air density calculation), 

assumed as 464° R 


With the use of these assumptions and arbitrary constants, equa- 
tions (l) to (7) may be combined as follows: 


Q. _ + +• 

ha " t s ‘ *0 - 
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The term q/ha of equation (9) is measured dimensionally in 
degrees of temperature, and is now expreseed as the arithmetic sum of 
five independent terms, as follows: 


5 - T 1 - T 2 + K ( T 3 

where from inspection of equation (9), 


T 1 = (*8 - ^oH 1 


V 

5X10- 
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/ e s\ 
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- T 4> + T 5 

+ MwOy ^ 




(10) 

(ID 

( 12 ) 

(13) 

(14) 

(15) 


These T terms are in a form that permits graphical presentation, as 
will he shown in the section Graphical Solution. 

Mass transfer . - The rate of evaporation from a wetted surface 
above 32° F is given in reference 2 as 


M ev = 0.622 K 


e s ^ e i 

\ p * c py 


(16) 


When equation (6) is substituted in the preceding equation, 


Mgy = 0.622 


V /e s fo N 
c p \i ~ p oy 


(17) 


or 


h K 

M ev - 1§66 < T 3 - T 4> 


( 18 ) 


When the surface is dry (no impingement or run back) , K and 
are zero and equation (10) simplifies to 
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A - Ti - T 2 + T 5 (19) 

When the rate of impingement My is lower than the rate of evapor- 
given by equation (18), the surface remains essentially dry 
because the water evaporates as it impinges, and the following equation 
applies : 

£ = T 1 ~ r 2 + t 5 + 1066 ^ ( 20 ) 

Equation (20) results from the addition to equation (19) of a term 
which includes the heat required to evaporate the impinging water 
completely. 


Graphical Solution 

General solutions . - The general equations of heat and mass trans- 
fer from a surface' subjected to water impingement and heated above 
freezing (equations (10), (18), and (20)) have been formulated using the 
five terms T]_ to T 5 . These T terms are presented graphically in 
figure 1 and are determined by intersections, as follows: 


Abscissa 


Sloping line 

T]_ = Intersection 


and 

"^s ~ t 0 

Tg = Intersection 

“vAi 

and 

O 

> 

t 3 = Intersection 

Pz 

and 

^s 

T4 = Intersection 

Po 

and 

t 0 

Tg = Intersection 

P Z 

and 

PO 


A numerical example is shown by arrows and dashed lines on figure 1 
to illustrate the use of the graph for the two principal solutions. The 
calculations required for the two solutions (using equations (lo) and 
(18)) are quite simple after the eight initial variables are given. 

Figure 1 may be used in conjunction with equations (10) and (18) 
to solve for any of the factors appearing therein, provided only one 
factor is unknown. Solutions may be made directly for the factors q, 
M ev> K > and V Q , whereas the other factors, h a , p Q , p z , t g , and t Q each 

appear twice and must be obtained by trial- and-err or. These trial- 
and-error solutions converge rapidly because the calculations required 
are minor and trends can be visualized from the graphs in figure 1. 

When q is zero, a solution for t s yields the unheated datum,' or 
equilibrium, surface temperature. 
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The coordinate scales used in the working chart of figure 1 should 
he large enough for the usual requirements of accuracy, permitting anti- 
icing calculations to be made quickly from the graph alone. In addition, 
the graph illustrates the relative magnitude of the various component 
terms and facilitates comparisons, check calculations, and evaluation of 
assumptions and short-cuts. 

The ambient temperature, pressure, and velocity terms in equa- 
tion (10) are deter min ed for the undisturbed stream conditions just 
ahead of the body in question. When the body is located in a pipe or 
duct, the ambient terms should be evaluated for conditions in the duct 
flhpari of the body, and not for the external air -stream conditions. 

Point and area solutions. - The solutions described previously have 
been based on local point values of the primary variables and expressed 
for unit surface areas. If appreciable variations occur in the surface- 
wise distribution of these primary factors, the local q and Me V solu- 
tions, as previously determined, should be obtained for several points 
along the surface in the direction of flow starting at the stagnation 
region. The total heat and mass transfer from a specific surface area 
may then be obtained by graphical or numerical integration. In this 
point -by-point manner, abrupt changes in surface conditions will be taken 
into consideration, for example, (a) the decrease in the external heat- 
transfer coefficient immediately downstream of the stagnation region, 

(b) the downstream limit of water impingement, (c) the transition from 
i fl TnirflT to turbulent heat -transfer, (d) the aft limit of run-back rivu- 
lets, and (e) any sudden changes in surface contour or surface tempera- 
ture. Furthermore, in this manner it is possible to determine how far 
downstream on a surface the impinged water will run back before it is 
evaporated by summing up the cumulative impingement and evaporation from 
the stagnation region downstream to a particular point; thus 


n n 

^ ~ M ev ) ££ 

S=0 


( 21 ) 


When Mj. becomes zero all the run back will have been evaporated. 

Although the point -by-point solutions discussed are usually required 
for accuracy, the graph and associated equations are equally applicable 
to area solutions in which the primary factors are constant over a given 
area, or approximately constant and can be conveniently averaged'. In 
such a case, area solutions are obtained in one operation by multiplying 
q and Me V by the specific areas involved. 

Solutions with gas heating. - With a heating source of hot gas, the 
local heating rate q is not as directly measurable as with electric 
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heating; whereas the input wattage is known for an electric heater, with 
gas heating, gas temperature , flow rate, and passage design are usually- 
known. For thin conductive skins between the gas supply and the external 
air stream and neglecting radiation, the following equation may be used: 

<1 = h g (t g - t s ) (22) 

The q values of equations (22) and (10) may then be equated and solved 
for any one unknown. If, for example, surface temperature t s is 
unknown, a trial solution may be made by estimating t s and proceeding 
until a check is obtained on the q values of equations (22) and (10). 

Usually an appreciable drop in gas temperature tg occurs in the 
direction of flow, necessitating a stepwise analysis in the flow direc- 
tion. Detailed treatments of surface evaporation and gas-heating 
considerations are given in references 1 and 4 . 


SUPPLEMENTARY GRAPHS 

Of the eight factors and variables required for. solution of equa- 
tions (10) and (18) (with fig. l) , several are involved and tedious to 
calculate. Graphs that simplify these calculations have been included 
wherever possible to do so in general terms. In the usual icing or 
heat-transfer problem, the ambient conditions, the free-stream velocity, 
and the geometric-shape factors are known and constant for a given set 
of calculations; these known factors can sometimes be plotted to yield 
the required factors, which may be considered dependent variables. 


Local Pressure and Velocity 

The conversion between local static pressure pj and local stream 
velocity vj for the subsonic range of Mach number M is shown in fig- 
ure 2 for compressible dry-air flow. Local velocity ratios may be 
directly converted into local pressure ratios after the free-stream Mach 
number curve has been located. The free-stream Mach number is determined 
from the relation 

V 0 

Mo = (23) 

49.1 //Tq 


Heat-Transfer Coefficients. 

External air flow . - In general, the value of the dry-air convective 
heat-transfer coefficient h a is difficult to establish and several 
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important assumptions are necessary. Heat-transfer coefficients are 
available for air flow over specific shapes, such as flat plates and 
cylinders (reference 6), over airfoils (references 1, 3, 7, and 8), and 
over bodies of revolution (references 9 and 10) . The coefficients for 
flat plates and cylinders are, at present, general in application and 
are presented graphically in figure 3. In figure '3(a) the heat-transfer 
coefficients are given for laminar and turbulent air flow over flat 
plates and for the leading-edge stagnation point of cylinders. The 
ratio of local to stagnation point heat-transfer coefficient over the 
forward surface of a cylinder is shown in figure 3(b) . The coefficients 
(fig. 3(a)) are determined by the use of an air -flow parameter consist- 
ing of air velocity, static pressure, and a significant length dimension. 
In order to simplify these relations, the absolute -temperature terms for 
surface and ambient air were taken as the same constants assumed for the 
derivation of equation (10). Errors resulting from this simplification 
are negligible when the many uncertain factors present in icing condi- 
tions are considered. 

The transition from laminar to turbulent heat -transfer is not 
easily located when the surface is covered with a water film or rivulets* 
In such cases, transition is very apt to be forward of the transition 
point- located for flow over a dry surface. Limited evidence (refer- 
ence 3) indicates that transition probably occurs near the downstream 
limit of droplet impingement on an airfoil. For unusual shapes and con- 
ditions, estimates are required for locating transition. 

Internal gas flow. - The heat-transfer coefficient between a gas 
flowing in a confined passage and the wall of the passage is a complex 
function. The coefficient is known to vary with the distance from the 
entrance to the passage, with the passage shape, with position around 
the passage perimeter, with the gas -flow rate, and with the gas tempera- 
ture (references 6, 11, and 12). The flow condition generally applicable 
to anti-icing calculations is that of fully developed turbulent flow in 
passages for which entrance length and position along the perimeter are 
assumed to have no effect. A condition of this type is shown in fig- 
ure 4, in which the internal heat-transfer coefficient is obtained from 
the parameter hg ? ^ (Ap/p0‘ 2) plotted as a function of the gas flow and 
gas temperature. The term (Ap/P°‘ 2 ) varies only with passage geometry 
and will generally be constant for a series of calculations involving a 
specific flow passage. Figure 4 does not include radiant heat transfer. 


Water -Interception Rate 

Sufficient data are not available to permit generalization of the 
impingement term My- for application to airfoils and bodies at various 
flight and ambient conditions. Impingement values therefore must be 
estimated on the basis of published data for specific airfoils (refer- 
ences 13 and 14) and for cylinders and spheres (reference 15) and from 
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the knowledge of impingement limits learned from icing observations and 
photographs (references 3 and 16 to 20). Occasionally, assistance is 
obtained from the maximum condition of impingement that occurs with high 
airspeed, large droplet size, and small body size, for which straight 
line trajectories of water drops may be assumed. Under these conditions 
Mw = (m/4.45 ) (Vq sin t), where T is the angle between the tangent to 
the surface at the point of impingement and the free-stream relative 
wind vector. 


Surface Wetness 

In the absence of exact data, the wetness fraction K in the basic 
equations is usually taken as 1 in the impingement area. Near the 
downstream limit of impingement, rivulets of run-back water commence 
and K decreases very rapidly. The K factors in regions of rivulet 
flow probably vary with many conditions, but limited data (reference 3) 
indicate that the value of K decreases in this region from a value 
near 0-3 down toward zero at the point of complete evaporation of run 
back. 


CONCLUDING REMARKS 

In the preceding derivations, the assumption is made that the 
moisture in the ambient air does not have time to change state during 
the adiabatic pressure and temperature changes in the local air stream. 
Because of the high speeds and short distances involved in most anti- 
icing calculations, this assumption appears reasonable. Answers obtained 
by use of this assumption do not differ appreciably from those using the 
assumption of phase changes in equilibrium with the local air tempera- 
ture. The difference between the two solutions is greatest when P 7 /P 0 
is most removed from unity and when low heating rates are used. 

Errors in the present graphical solution may become appreciable at 
high subsonic speeds because an incompressible-flow relation was used in 
eliminating the velocity-ratio term from the equations. Recovery factors 
for wetted surfaces at high speeds may also deviate from the value 
assumed herein. The errors in solution due to the temperature factors 
that were taken as constants near the mean of their expected ranges are 
in the order of 3 percent for the range of graph presented. The fore- 
going errors are negligible, however, compared with the uncertainty and 
the latitude involved in estimating values of the basic variables, 
especially the water-interception rate and the heat-transfer coefficient. 
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The form of solution presented herein is intentionally made general 
for maximum usage. The solution is applicable to practically any body 
shape for which the flow conditions are known, and is not dependent upon 
any specific assumptions as to body geometry, location of transition, 
wetness of surface, and method or extent of heating. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 11, 1952 
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APPENDIX - SYMBOLS 


C P 

c w 

D 

e 

g 

h 

J 

K 

L 

M 

Mev 

M r 

m 

P 

P 

<1 

R 

S 

T 


The following symbols sire used in this report: 

cross-sectional area of duct or pipe passage, sq ft 
specific heat of gas at constant pressure, 0-24 Btu/(lb)(°F) 
specific heat of liquid water, 1.0 Btu/(lb)(°F) 
diameter of cylinder, ft 

partial pressure of water vapor, (corresponding to sat- 
urated air unless noted otherwise), in. Hg 

acceleration due to graylty, 32.2 ft/sec 2 

dry convective heat-transfer coefficient, Btu/(hr)(sq ft)(°F) 

mechanical equivalent of heat, 778 (ft)(lb)/Btu 
surf ace -wet ness fraction 

latent heat of vaporization of water, Btu/lb 
Mach number 

rate of evaporation of water, lb/(hr)(sq ft) 
rate of run back of surface water, lb/(hr)(ft span) 
rate of interception of water, lb/(hr)(sq ft) 
liquid -water content of ambient air, g/cu m 
inside perimeter of duct or pipe, ft 
absolute static pressure, in. Hg 
rate of heat transfer, Btu/(hr)(sq ft) 
kinetic energy recovery factor 

surface distance from stagnation point in stream direction, 
ft 

o 

absolute temperature, R 
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t 

V 

w 

w 

X 

r 

p 

T to T 
1 5 


temperature, °F 

relative velocity between surface and stream, ft/sec 

rate of heated gas flow, lb/sec 

rate of heated gas flow, lb/hr 

Hardy's evaporation factor (equation (3)) 

angle of impingement of water droplet upon surface, deg 

density, ( lb) (sec^)/ft^ 

temperature terms defined by equations (ll) to (15) and 
determined by intersections on figure 1, °F 




cylinder center -angle between stagnation point and point on 
cylinder surface, deg 


Subscripts: 


0 

a 

d 

g 

l 

lam 

s 

stag 

t 


undisturbed ambient stream 

external air 

heat -transfer datum 

internal hot gas 

local 

laminar flow 
surface 

cylinder stagnation point 
turbulent flow 


w 


water 
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(a) Dry-air heat- transfer coefficients for flat plates and cylinder stagnation points as functions of flow 

velocity, static pressure, and significant length. 



Ratio of local to stagnation point heat-transfer 
coefficients for cylinder, h a 7 /h„ 






Gas flow, W, Ib/sec 
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NACA - Langley Field, Va. 


Figure 4. - Determination of turbulent pipe-flow heat-transfer coefficient in terms of gas flow, temperature 

passage size. 





